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ABSTRACT 



We present a spectral analysis of the Far Ultraviolet Spectroscopic Explorer (FUSE) 
spectra of the eclipsing double-line spectroscopic binary EM Cygni, a Z Cam DN system. 
The spectrum, obtained in quiescence, consists of 4 individual exposures (orbits): 

two exposures, at orbital phases (p ^ 0.65 and (j) ~ 0.90, have a lower flux; and two 
exposures, at orbital phases (p = 0.15 and 0.45, have a relatively higher flux. The change 
of flux level as a function of the orbital phase is consistent with the stream material 
(flowing over and below the disk from the hot spot region to smaller radii) partially 
masking the white dwarf. We carry out a spectral analysis of the FUSE data, obtained at 
phase 0.45 (when the flux is maximal), using synthetic spectra generated with the codes 
TLUSTY and SYNSPEC. Using a single white dwarf spectral component, we obtain a 
white dwarf temperature of 40,000 ± lOOOK, rotating at lOOkm/s. The white dwarf, 
or conceivably, the material overflowing the disk rim, shows suprasolar abundances of 
silicon, sulphur and possibly nitrogen. Using a white dwarf+disk composite model, we 
obtain that the white dwarf temperature could be even as high as 50,000K, contributing 
more than 90% of the FUV flux, and the disk contributing less than 10% must have 
a mass accretion rate reaching lO^^'^M0/yr. The single white dwarf model fits the 
absorption lines better than the white dwarf+disk model, but the white dwarf+disk 
model fits better the continuum in the shorter wavelengths. In both cases, however, 
we obtain that the white dwarf temperature is much higher than previously estimated. 
We emphasize the importance of modeling the spectra of EM Cyg around phase < 



^Based on observations made with the NASA-CNES-CSA Far Uhraviolet Spectroscopic Explorer. FUSE is oper- 
ated for NASA by the Johns Hopkins University under NASA contract NAS5-32985 

^Visiting at the Space Telescope Science Institute, Baltimore, MD 21218; godon@stsci.edu 
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0.5, when the white dwarf and disk are facing the observer, and we suggest that the 
discrepancy between the present analysis and previous spectral analysis might be due 
to the occulting effect of the stream veiling the white dwarf and disk. 

Subject headings: accretion, accretion disks - novae, cataclysmic variables - stars: indi- 
vidual (EM Cyg) - ultraviolet: stars - white dwarfs 



1. Introduction 



1.1. Dwarf Novae 



Dwarf Novae (DNe) are a class of non-magnetic cataclysmic variables (CVs) in which a white 
dwarf (WD) accretes matter from a main-sequence star (the secondary) filling its Roche lobe by 
means of an accretion disk. DN systems are characterized mainly by their brightness variations, 
due to a periodic change in the accretion rate. Ongoing accretion at a low rate (quiescence) is 
interrupted every few weeks to months by intense accretion ( o utburst) of da ys to weeks which 
increases the luminosity of the system (iHack &: La Dous I Il993l : IWarner I Il995l ) . This outburst is 
believed to be due to an instability in the disk (the disk instability model - DIM (jCannizzo Ill998l )). 
and it is known as a dwarf nova accretion event (or dwarf nova outburst). 

Due to their long term accretion, the WDs in DNe have temperatures ~ 15, 000 — 50, OOOK 
(higher than field WDs), and therefore peak in the far-ultraviolet (FUV). The long periods of quies- 
cence allow for long exposure observations (~hrs) covering all the orbital phases (this is especially 
important for systems with a high inclination, as the WD can be eclipsed or masked at given orbital 
phases). Because of that, DNe, unlike other CVs, are ideal candidates to study the accreting WD, 
as during quiescence the WD is exposed in the FUV for long periods of time. FUV spectroscopic 
observations are used to determine the WD parameters such as temperature, gravity, rotational 
velocity and chemical abundances. The disk mass accretion rate of many systems can also be 
deduced accurately at given epochs of outburst or near outburst using spectral fitting techniques. 
In addition, DNe offer a fairly reliable e stimate of the i r distances via the ab solute magnitude at 
maximum versus orbital period relation (IWarner Ill995l : iHarrison et al. 112004 ). As a consequence, 
DNe are ideal candidates to observe, as one can more easily derive their system parameters, which 
in turn can be used to test the theories of CV evolution. 



Dwarf nova systems are further divided into sub-types: U Gem systems, SU UMa systems, and 
Z Cam systems. The U Gem systems are the typical DNe, i.e. with normal DN outbursts; the SU 
UMas exhibit both normal DN outbursts, and superoutbursts, which are both longer in duration 
and higher in luminosity than normal DN outbursts; an d the Z Cam systems have standstills where 
they remain in a semi-outburst state for a long time (IRitter &: Kolb II2003I ). The binary orbital 
period in CV systems ranges from a fraction of an hour to about a day; however, there is a gap in 
the orbital period between 2 and 3 hours where almost no CV systems are found: the period gap. 
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U Gem and Z Cam DN systems are found above the period gap, while the SU UMa DN systems 
are found below the period gap. 

In this paper we present the FUV spectral analysis of the eclipsing Z Cam system EM Cygni 
observed during quiescence with FUSE. 



1.2. EM Cygni System Parameters 



First suspected to be related to CVs bv iBurbidge Sz Burbidge I (|l953l ). EM Cy g is a bright, 

doubl e- lined spectroscopic binary variable ([Kraft I ll964l l. Its eclipses were discovered bv lMumford Krezminsk: . 
(|l969l ) and the orbital period was found to be almost 7h (0.29 90912 day). Its pe riod was expected 
to decrease (by ~ 10^^^ — 1 0~^^ days/cycle; Pringle ( 1975 1: Mumford ( 1980 ). however it has 



1984 : 



Csizmadia et al. 



2008) with an upper 



been found to be constant (jBeuermann &: Pakull 
limit P/P < 2.3 x 10~^^/days. Its outburst cycle is of the order of 25 days and rapid coherent 
oscillations are detected dur i ng th e late stages of outbursts with periods of the order of 14.6-16.5 
seconds (see IStiening et al. I ()l982l ) and references therein). 



EM Cyg is so far the only known Z Cam subtype eclipsing DN system. The eclipses are 
important as they allow, in theory at least, for a detailed study of the system. In addition, because 



studied (e. g 



Patterson (1981 



Robinson 



Stover et al. 



Nevo k Sadeh 


( 


1978) 


(1981): 


Szkodv 1 ( 


1981 


); 



the determination of the stellar masses of the system by iRobinson I ()197 



was effected by the presence of a third star, Arkadash]^, along the line of sight (North et al. 



Jameson et al. ( 1981 



see also 


Stienine et a 




( 


1982)). However, 


Robinson (19741) and 


Stover et al. 


1981) 



This third star is not related to the system. The most recent study (IWelsh et al. 
a mass ratio q = 0.7 7 and a WD mass M^^ = IMq. The distance found by 



200C 



2OO7I) indicates 



Bailev I (|l98ll ) 



Jameson et al. I (jl981), d ^320-35 pc, c ould also be affected by the contribution of Arkadash 
in the red light, and I Welsh et al. I (j2005l ) pointed out that d could be as large as ^450-500pc . 



However, using the rn aximum magnitude at outburst ve rsus orbital period re lation (jWarner 
Harrison et al. II2004I ). one finds a distance of ~ 400pc (I Winter &: Sion II2003I ). 



1995 



A recent broad-band photometric study of the system (jSpogli et al. I l2003l . l2005l ) shows the 
magnitude of EM Cyg varies between V=13.3 at minimum and V=12.5 at maximum, it reaches 
V=14.4 during the eclipse, and V=12.9 in standstill (note, however, that data from the AAVSO 
indicates it can reach V~12.0). EM Cyg is located at a very low Galactic latitude (6 = 4.28°), 
where the Galactic reddening reaches almost 0.5. However, the reddening towards EM Cyg is only 



^ The third star was named Arkada§ ("friend" in Turkish) bv IWelsh et al. I ( 20051 ) . in memory of Janet (Han- 
nula) Akyiiz Mattei (1943-2004). Janet Mattel determined an outburs t cycle of 25 days for EM CyR, and p r ovide d 
AAVSO data and information on the state of EM Cyg many times (e.g. Stover et al. ( 198lh : Stiening et al. ( 1982h : 
North et al. I (|200Q )). as she did for many variable stars on countless occasions. 
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0.05 ±0.03 (Verbunt 



1987 



La Dous 



199ll : iBruch &: Engel 1 1 19941 ) . The system has an inchnation 



i of only 67 it 2deg (jNorth et al. Il200d ). however, because of its large disk and secondary star, the 
disk is partially eclipsed, while the WD itself is not eclipsed. The geometrical configuration of the 
disk, WD and secondary of the system is shown in Figure 1 at phase (p = 0/1.0. All the system 
parameters are recapitulated in Table 1. 



1.3. Previous Spectral Analysis of EM Cyg 



EM Cyg was observed in quiescence with FUSE by I Welsh et al. I (|2005l ). who carried out a 
spectral fit after combining the FUSE spectrum with an lUE spectrum. Their spectral fit consists 
of a disk with M = 2.5 x W'^'^Mq/yv and a O.SSM© WD with T w 10 - 20, OOOK, assuming 
a distance of 350pc. The model did not fit the flux in the short wavelengths of FUSE which is 
indicative of a much higher temperature. However at T > 20, OOOK the model spectrum became 
too blue and the fit was noticeably worse in the lUE range. 



Winter fc Sion I ([2003) carried out a spectral fit to the lUE spectrum SWP 08088 of EM Cyg 
taken in quiescence (with a flux of ~ 3 x 10~^'*ergs cm~'^A~^) and obtained that the best fit is a 
disk with M = 5 x lO^^^/o/yr and a WD of IM© wit h T < 24, OOOK (cont ributing only 8% of 
the flux) with a distance of 210pc and i = 75°. Later on, lUrban &: Sion I (|2006l ) analyzed the same 
lUE spectrum and obtained the same results except a d istance of 350pc, when dereddening the 
spectrum assuming E(B-V)=0.05. [Hamilton et al. I (120071 ) analyzed the /[/£' spectrum SWP 07297 
of EM Cyg in outburst assuming d=350pc, i = 75°, and obtained that the flux comes mostly from 
a disk with M = 7 x 10"^ M© /year. 



In the spectral analysis of IWinter &: Sion I (j2003l ) and lUrban &: Sion I (120061) ori l y the lUE 
spectra of EM Cyg in quiescence was modeled. While in the work of I Welsh et al. I (j2005l ) the 
fitting of the s hort wavelength rang e of FUSE (< lOOOA) was not carried out satisfactorily (as 
pointed out by I Welsh et al. I (|2005l ) themselves). It is the purpose of this work to carry out a 
spectral fit to the FUSE spectrum of EM Cyg in quiescence including the short wavelength range, 
and to obtain the physical parameters of the exposed WD. 

We give particular importance to the orbital phase at which the individual FUSE exposures 
were obtained. We show that at particular orbital phases the WD might be veiled by the stream 
from the first Lagrangian point overfiowing the disk. We remark that previous estimates might have 
modeled a flux attenuated by as much as 75%, and consequently obtained a low WD temperature. 
We model here the FUSE spectrum of the WD when veiling is minimum and obtain a temperature 
T ~ 40, OOOK and possibly as high as 50,000K. 

In the next section we present the four individual FUSE exposures and compare them with 
the lUE spectra taken in quiescence; in section 3 we review the stream-disk interaction and in 
particular we check the possibility of the stream overflowing the disk; in section 4 we describe the 
synthetic stellar spectral code as well as our fltting technique; the results are presented in section 
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5; and a discussion with our conclusions are given in the last section. 



2. Observations of EM Cyg 

2.1. The FUSE Spectrum 

The FUSE data of EM Cyg were obtained on 2002 Sep 5 (HJD2452522-3, when the system 
was in a low state, see Figure 2) through the 30"x30" LWRS Large Square Aperture i n TIME TAG 



mode. The data were processed with the l atest and final v ersion of CalFUSE (v3.2.0; iDixon et al 



(|2007l )). We process the FUSE data as in iGodon et al. I ([200a), and details of the procedure can 



be found there. 

The FUSE spectrum of EM Cyg consists of 4 individual exposures (one per orbit). For the 
purpose of line identifications we present the averaged spectrum in Figure 3. We identify many ISM 
molecular hydrogen absorption lines as well as some of the metal lines usually seen in the FUSE 
spectra of accreting WDs. The metal lines we identify are as follows: 

- the Sili (1021 A), Siv (1062 k 1073 A), S vi (945 A) lines and Ovi doublet, all indicative of a 
higher temperature T ~ 35, OOOK and much higher; 

- the Sim (1114 A ), Siiv (1122 A ), and Nil (1085 A ) lines, sometimes indicative of a lower 
temperature T ~ 25, OOOK, but not inconsistent with a higher temperature; 

- the Cm (1175 A) absorption line, present at all WD temperatures, as long as the carbon abun- 
dance is a fraction of its solar value. 

The C 11 (IOIOA& 1066A) absorption lines, seen at lower temperatures (T ~ 25 — 30, OOOK), are not 
detected here (though the Cii (lOlOA) could be contaminated with ISM molecular absorption). 
There are some broad emission features from C iii (977 A), and the O vi doublet. The C iii (1175 A) 
presents a P-Cygni profile. The sharper absorption lines are from the ISM and possibly also from 
circumbinary material. It is likely that the nitrogen absorption lines (N i & N li) are contaminated 
with terrestrial absorption. All these lines are marked and identified in Figure 3. There are a 
few additional absorption features/lines that we cannot identify, mainly in the longer wavelengths 
(these lines are not labeled in the lower panel of Figure 3). The sharp emission lines are from air 
glow (geo- and helio-coronal in origin). The broad Ly/3 absorption feature is not clearly seen, point- 
ing to a high temperature, some broad oxygen doublet emission contamination, and/or possibly 
a rotational velocity broadening from a disk component. In the short wavelengths (< 950A) the 
flux does not go to zero and the continuum there is about 1/3 of the flux at longer wavelengths 
(> lOOOA), indicative of a high temperature and consistent with the presence of higher ionization 
metal lines. 

As mentioned, this FUSE spectrum consists of 4 individual exposures taken during 4 consec- 
utive orbits of the FUSE telescope. The four exposures corr espond to the binary pha ses ~0.45, 



0.65, 0.90, and 0.15, where (/> = corresponds to the eclipse (jWelsh et al. 1120051 . 120071 ). The obser- 



vations log for the individual FUSE exposures is given in Table 2. Using an analysis of the lines 
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in the four individual FUSE exposures, IWelsh et al. I (|2005l . 120071 ) explicitly showed how the radial 



velocity changes as a function of the orbital phase. However, they did not consider the variation of 
the continuum flux level from exposure to exposure, i.e. the variation of the FUV flux as a function 
of the orbital phase. In Figure 4 we show how the flux varies from exposure to exposure, and this 
change is also registered in the visual magnitude of the system as shown in data from the AAVSO 
in Figure 2. We find that at phases 0.15 & 0.45 (exposures 4 & 1) the FUSE continuum flux is 
larger (by up to a factor three) than at phases 0.65 & 0.90 (exposures 2 & 3). 

In order to understand the variation of the FUV flux, in Figures 5 we show the binary config- 
uration at the times the 4 FUSE exposures were obtained. From the orbital configuration we see 
that the flux in exposures 2 & 3 is attenuated when the hot spot region and the stream overflow 
region are in front of the observer. In the next subsection we carry out an examination of the lUE 
spectra of EM Cyg, and we find that they too present some flux variation as a function of the 
orbital phase. We suspect that the FUV flux is possibly attenuated by absorption from the stream 
material (originating at LI) overflowing the disk edge and reaching a height above (and below) the 
disk of the order of z/r ~ 0.4 We further check this suggestion in section 3. 



2.2. The lUE Spectra 



Since the FUSE Qux varies as a function of the orbital phase it is interesting to compare it with 
other spectra obtained in quiescence; this will also help us interpret results from previous spectral 
analysis. The existing archival HST/STIS spectra of EM Cyg were taken during outburst, just a 
week before the acquisition of the quiescent FUSE spectrum (see Figure 2), and the three existing 
i^f/T spectra were also taken in outburst. Because of that we disregard the STIS and HUT spectra. 



However, some of the e xisting lUE spectra were taken in c 



uiescence, and later were modeled by 
([200a). We present and discuss these 



Winter fc Sion I (120031 ^: IWelsh et al. I ([20051); lUrban fc Sion 
spectra here. 

In Figure 6 we compare three lUE SWF spectra taken in quiescence together with the FUSE 
exposures in the spectral range ~ 900 - 1,400 A . Since FUSE exposures 2 and 3 have about the 
same flux, we combine them together, and the same for exposures 1 and 4. We see that SWF 08088 
(taken in 1980) has a much lower flux, similar to exposures 2 and 3. SWP 17591 (taken in 1982) 
has a flux level comparable with exposures 1 and 4, while SWP 17592 (taken 90min after SWP 
17591) has a higher flux level. 

The typical quiescent V of the system is 13.5 bu t it has been observed as faint as 14.5 in 
quiescence (possibly at eclipse (ISpogli et al. Il2003l . l2005l )). From validated AAVSO data, it appears 
that, at the time of the FUSE observations, the system had a magnitude of ^ ~ 13.1 — 13.9 
(Figure 2), namely EM Cyg wasn't in its lowest possible quiescent state at the time of the FUSE 
observations. This is clearly seen in Figure 2, where the previous quiescent state on JD2452505 is 
about 0.5mag fainter than the quiescent state during which the FUSE observations (JD2452522) 
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were carried out. The lUE SWP 17591 and SWP 17592 spectra were obtained during a brighter 
state, when EM Cyg was on the rise to outburst, its magnitude around 13.2-13.5 (but still lower 
than a standst ill with V ~ 12.5 - 13.5). The lUE SWP 08088 spectrum was at F ~ 13.6 - 13.7 



(|Szkodv 111981 



The three lUE SWP spectra are complemented with three LWR spectra. All the quiescent 
(LWR and SWP) /f/i? spectra are listed in Table 2 together with their corresponding orbital phase. 
The increase in flux between SWP 17591 and SWP 17592 over a time scale of 90min is most 
probably due to orbital variation (as in FUSE), and the matching LWR spectra 13862 and 13863 
do not show any variation. As to the SWP 08088 and (matching) LWR 07056 spectra, their lower 
flux is due to a deeper quiescence and to orbital variation too. 

Because of these changes in flux level (due to both orbital variation and intrinsic changes in 
the sour ce) combining the F USE and lUE spectra together is tricky and this might explain the 



difficultv IWelsh et al. I (|2005l ) were having to match both the short wavelengths of FUSE and the 
longer wavelengths of lUE'va. the same flt. As a consequence we do not combine the F[/S'£' spectrum 
of EM Cyg with any lUE spectrum. The shorter wavelength range as well as the higher S/N of 
FUSE over lUE makes the FUSE spectrum of EM Cyg ideal for determining the temperature of 
the accreting WD. 

From the FUSE & lUE spectra of EM Cyg, we found that the FUV flux between ~900A and 
2000A is attenuated around phase (f) ~ 0.6 — 0.9. Since the FUSE spectrum presents many deep 
absorption lines, it is tempting to suggest the possibility of an absorbing component moving in front 
of the observer around phase ~ 0.75, as the binary rotates. We further investigate this possibility 
in the next section. 



The Stream-Disk Interaction 



The only components of the system that could be in the fleld of view of the WD/inner disk 
around (p ~ 0.6 — 0.9 are the region downstream the hot spot and the stream possibly overflowing 
disk at the hot spot. In order to check this possibility, in this section, we consider previous studies 
of the stream-disk interaction, and also present a simple ballistic trajectory simulation using the 
binary parameters of EM Cyg. We show that for the binary parameters of EM Cyg, stream material 
could be found high enough above the disk of the plane, around phase (j) ~ 0.6 — 0.9 to actually 
veil the WD and inner disk. 



Pioneering analytical studies by lLubow &: Shu I (jl976l ): lLubow I (jl989l ) have shown how part of 
the stream material possibly flows over and under the disk and falls back onto the disk at smaller 
radii at an angle 4>impact ~ 

140 - 150° ((b 

impact is measured counter clockwise from the line joining 
the two stars) , corresponding to an orbital phase (j) ~ 0.6. Observationally, Doppler tomography 
(jMarsh Ill985l ) has been used for a quarter of a century to conf irm the s treain disk interaction and 
the impact at smaller radii in some system (e.g. Z Cha). In iLubow I (jl989l ) some of the stream 
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matter was flowing over the disk at a scale height ~ O.Ola, where a is the binary separation. For 
a disk of diameter Rdisk = 0.3a, this corresponds to a scale height of z/r 0.03 at the edge of the 
disk. Namely, in these analytical studies, it was obtained that stream material overflows the disk 
rim (at the hot spot) and continue toward smaller radii while basically grazing the disk's surfaces 
(at a few vertical scale lengths - a few ^ H). Therefore, based on these studies alone, in order 
for the material to veil the WD, the inclination of the system has to be cosi ~ 0.03 ^ i ?a 88° or 
larger: the system has to be nearly edge-on. 



However, more recent full 3D hydrodynamical simulations, e.g. ([Blondin 1119981 : iKunze et al 



200ll : iBisikalo et al. II2003I ). have followed in detail the interaction of the material streaming from 
the first Lagrangian point (LI) with the disk. The main result from these simulations is that the 
stream does not simply grazed the disk's surfaces, but matter is actually deflected vertically from 
the hot spot and flows in a more or less diffuse stream to inner parts of the disk, hitting the disk 
surface close to the circularization radius rc (see iLubow Sz Shu I (jl975l ) for details on Vc) at orbital 
phase 0.5. In these simulations, the stream material deflected vertically at the hot spot region 
can reach an altitude z ~ r. This stream deflection is believed to cau se X-ray absorption in CVs 



(and LMXBs) around orbital phase 0.7, if the inclination is at least 65° (jKunze et al. Il200ll ). These 



simulations also explain stre ams of material p e rpend icular to the orbital plane as observed in Algol- 
type binaries (e.g. U CrB; lAgafonov et al. I (|2009l ) and references therein). These results imply 
that stream material can be deflected vertically, and reaches a height z/r = cosi = cos 65° = 0.42 
and lands onto the disk at (pimpact = 180° (corresponding to an orbital phase 0.5) (jKunze et al. 



2001 



We suggest the same might be happening in EM Cyg, and the decreased flux around (j) ~ 0.7 
might be caused by the deflected stream flow. In order to verify the plausibility of this scenario, 
we simulate the stream flow from the LI point using the ballistic trajectory approximation, and we 
deflect the ballistic particle at the hot spot by an angle 6 (measure from the plane of the disk). Since 
the stream flow consists of a infinite number (in the continuity limit) of particles (and since the 
stream has a finite thickness extending vertically), the ballistic trajectory is solved for all possible 
value of 6 (from —90° to +90°). A maximum altitude is reached for a critical angle of deflection 
9 = Ocri- This maximum altitudes depends on the values of the system parameters: veiling will not 
occur for every mass ratio, disk radius and inclination (the full de tails of the sim ulations as well as 
results for the entire space parameters of CVs are given elsewhere (jGodon II2009I )). Figure 8 depicts 
such a simulation for EM Cyg, where we show with arrows the 4 Ff/S'E' exposures. The mass ratio is 
0.77 and Rdisk = 0.3a. The stream is shown in red and the blue area represents that portion of the 
trajectory for which the material reaches z/r > 0.4, corresponding to veiling at i = 67°. Here the 
stream is deflected at the hot spot with an angle 9cri ~ 55° above (and below) the plane of the disk. 
This is the angle for which the altitude of the ballistic particle reaches a maximum. For a smaller, 
or larger, deflection angle 9, the blue region shrinks. For a smaller disk radius {Rdisk < 0.3) 
the maximum elevation above the disk increases, however, the hot spot moves counterclockwise 
towards (p ~ 0.75 and as a consequence exposure (3) is not as veiled as for Rdisk = 0.3a. This 
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points to the importance of the large disk in EM Cyg for deflecting material off the disk plane 
around 4) ~ 0.6 — 0.9. This picture is consistent with the full 3D hydrodynamical simulations of 

For different system parameters (inclination, mass ratio and disk radius) 



Kunze et al. 



(|200ll 1. 

the deflected material does not especially veil the WD, as the maximum vertical elevation of the 
ballistic particle does not always reach z/r = cosi. 

From this specific ballistic trajectory simulation with the system parameters of EM Cyg, we 
conclude that the lower flux in the FUSE exposures 2 and 3 can be accounted for with stream 
material being deflected vertically at the hot spot forming a veil in front of the WD between orbital 
phases ~ 0.65 and ~ 0.90 (for an inclination of 67°). The increase in flux between the /f/i? spectra 
SWF 17591 {(j) ~ 0.75) and SWF 17592 {(j) ~ 0.97) over a time scale of 90min can also be explained 
by veiling of the WD/inner disk by the stream overflow around phase 0.75 (keeping in mind that 
the WD/inner disk are not eclipsed around ~ 1.0). 

Furthermore, the stream material overflowing the disk edge lands onto the disk at a smaller 
radius correspondin g to orbital phase of up to ~ 0.5, and produces a bright spot there ([Marsh 



19851 : iLubow Ill989l ). Together with the hot spot, this contributes to additional flux in the longer 
wavelengths (LWR spectra). The emission from the hot spot and point of impact is possibly seen 
in the LWR spectra (13862, 13863) shown in Figure 7 as a plateau in the longer wavelengths 
(A ~ 2250 - 3000A), where one rather expects the flux to decrease ther e. This f l ux in the longer 
wavelengths was already associated with emission from the hot spot by ISzkody I (Il98ll ). Though 
we have no way to confirm that the plateau in the longer wavelengths of lUE is due in part to the 
impact point of the stream overflow, it is, however, consistent with this possibility. The plateau 
seen in the LWR spectra, is also contributing significantly to the longer wavelengths of the SWF 
spectra. This is especially true for the SWF 08088 spectrum. These are additional reasons not to 
model the lUE spectra, nor to combine them with the existing FUSE spectrum, as the hot spot 
and point of impact contribute flux to the longer wavelengths of the lUE SWF spectra. 



4. Synthetic Spectral Modeling 
4.1. Modeling the Stellar Spectrum 

We crea te model spectra for high-gravity stel lar atmospheres using the codes TLUSTY and 



SYNSFEcl (jHubeny lll988l : lHubeny Lanz Ill995l ). We generate photospheric models for a LOM© 



WD with effective temperatures ranging from 12,000K to 75,000K in increments of about 10 percent 
(e.g. 1,000K for T«15,000K and 5,000K for Tf«70,000K). We vary the stellar rotational velocity 
Vrotsin{i) from 50km s""*^ to 500km s""*^ in steps of 10 — 50km s""*^. 



We also use disk models from the grid of models generated by I Wade &: Hubeny I (jl998l ). We 



^ |http://nova.astro. umd.edu 



TLUSTY version 200, SYNSPEC version 48 
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select their models with a WD mass M = IMq and inclinations i = 60° & 75°. These models 
have solar abundances, but because of the velocity broadening the absorption lines are completely 
smoothed out, such that the disk models (i ~ 60° — 75°) are not sensitive to the chemical abun- 
dances. These disk models are either used as such or in combination with WD models we generate. 



After having generated grids of models, we use FIT ( Press et al.lll992l ). a minimization 
routine, to compute the reduced xt ix^ P^^ number of degrees of freedom) and scale factor values 
for each model fit. Before carrying out each synthetic spectral fit of the spectrum, we masked 
portions of the spectrum with strong emission lines, and air glow. While we use a minimization 
technique, we do not blindly select the least models, but we also examine the models that best 
fit some of the features such as absorption lines. 

The WD rotation (V^ofsm(i)) rate is determined by fitting the WD model to the spectrum 
while paying careful attention to the line profiles. We did not carry out separate fits to individual 
lines but rather tried to fit the lines and continuum in the same fit. We first fit solar abundance 
models and only then we change the abundances of the species. 



4.2. Modeling the ISM Hydrogen Absorption Lines 



We identify ISM lines in Figure 3 to avoid confusing them with the WD lines. The ISM lines 
are deep and broad and we decided to model them, especially since some of the WD lines (such as 
S IV AA1062.6 & 1073) are located at almost the same wavelengths. 

We model the ISM hydrogen absorption lines to assess the atomic and molecular column 
densities. This enables us to improve the WD spectral fit. The ISM models (transition tables) are 
generated using a program developed by P.E. Barrett. This program uses a custom spectral fitting 
package to estimate the temperature and density of the interstellar absorption lines of atomic and 
molecular hydrogen. The ISM model assumes that the temperature, bulk velocity, and turbulent 
velocity of the medium are the same for all atomic and molecular species, whereas the densities 
of atomic and molecular hydrogen, and the ratios of deuterium to hydrogen and m etals (inc l uding 
heliu m) to hydrogen can be a djusted independently . The model uses atomic data of iMorton I (|200d . 
20031 ) and molecular data of lAbgrall et al. I (120001). The optica l depth calculations of molecular 
hydrogen have been checked against those of iMcCandliss I (|2003l ). 



The ratios of metals to hydrogen and deuterium to hydrogen are fixed at and 2 x 10^^, 
respectively, because of the low signal-to-noise ratio data. The wings of the atomic lines are 
used to estimate the density of atomic hydrogen and the depth of the unsaturated molecular lines 
for molecular hydrogen. The temperature and turbulent velocity of the medium are primarily 
determined from the lines of molecular hydrogen when the ISM temperatures are < 250K. 

The ISM absorption features are best modeled and displayed when the theoretical ISM model 
(transmission values) is combined with a synthetic spectrum for the object (namely a WD synthetic 
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spectrum or a WD+disk composite spectrum). 



5. Results 



Though the distance es t imate to EM Cyg is most hke around 400pc ( Bailey lll98ll : IJameson et al 



1981 



Winter &: Sion I l2003l : IWelsh et al. 1120051^. w e do not fix the distance in our model fitting. 



We adopt a WD mass of IM© (jWelsh et al. 1120071 ) in all our models, with a corresponding radius 



of ~ 6,000km (giying Log{g) ~ 8.6). All the model fits are carried out on the first exposure (with 
the maximum flux). The reddening toward EM Cyg is 0.05 it 0.03. As a consequence, we model 
EM Cyg assuming both E(B-V)=0.0 and E(B-V)=0.05. We first carry out the fitting without 
dereddening the spectrum, and only later do we deredden the spectrum (assuming (E(B-V)=0.05) 
to assess the effect of reddening on the results. 

We first run single disk models. We find that in order to fit the short wavelengths of FUSE, the 
best single disk model has a mass accretion rate M ~ lO^^MQyr"^, i.e. far larger than expected 
for quiescence, and i ~ 60° — 75°. The i = 75° model is shown in Figure 9. It gives a distance of 
441pc and xl = 0.6916 (ah our X values are smaller than one due to the low S/N of the FUSE 
exposures). On top of the model, we add an ISM model in order to optimize the fit. The ISM has 
the following parameters Nh = 2 x 10^^ atoms/cm^ , = 1 x 10^'' molecules/cm^, T=160K, 
v=25km/s, D/H = 1 x 10^^, and Z=Ox solar. We did not try to optimize the ISM model itself, 
nor to obtain accurate parameter values for the ISM. The modeling the ISM is carried only in the 
context of improving the disk/WD model fit and to help identify the lines and spectral features. It 
is clear from Figure 9 that this disk model only vaguely follows the continuum of the spectrum and 
does not fit the absorption features (except for the ISM lines) . The two best single disk models are 
listed in Table 3. 

Next, we run single WD models. We find that in order to fit the continuum the WD must have 
a temperature of 40,000K, with a corresponding distance of 280pc . The x^ obtained is Xu ~ 0.58. 
The two first WD models in Table 3 shows how the distance and x^ change as a function of the 
WD temperature. In this first set we chose a guessed rotational velocity of 200km/s. In the next 
model (Table 3), we improve the x"^ value by varying the rotational velocity. The optimal value is 
somewhere between 90 and lOOkm/s, and we decide to fix the rotational velocity to lOOkm/s. 

In the next four models (Table 3) we change the abundances of silicon, sulphur and nitrogen 
while keeping all the other abundances to their solar value (we found that a solar carbon abundance 
is the best fit, though the main carbon line at 1175A has also some broad emission and makes it 
difficult to model). The best fit is obtained for enhanced abundances as follows: Sulphur: 10 
times solar. Nitrogen: 10 time solar and Silicon: 20-30 times solar. We obtain a distance less than 
anticipated: ~260pc. We present this best fit model in Figure 10. In order to optimize the fit we 
shift the synthetic WD spectrum to the blue by 0.5A , to account for the radial velocity shift. This 
improves slightly the x^ values (the values listed in Table 3 refiect this improvement). We also add 
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an ISM model to the WD model, however, the ISM model is not shifted. 

In the next model in Table 3 we deredden the spectrum assuming E(B-V)=0.05, and obtain 
the same results, except for the distance. Here the distance is barely 200pc. We note that assum- 
ing a WD mass of ~ 0.8Mq can increase the distance to 350pc without changing the remaining 
parameters of the fit. On the other side, adding an accretion disk does increase the distance too, 
and this is what we do next. 

For the WD+disk composite model, we first take our best single WD model and add a disk 
model to it, to try and improve the fit and increase the distance. We find that the mass accretion 
rate has to be at least M = lO~^'^M0yr~^ for the disk to contribute a noticeable flux. As the disk 
is added, the model deteriorates and give a larger Xu- We present such a model in Figure 11. The 
fit in the lower wavelength is not as good as the best fit single WD model presented in Figure 10. 

In order to try and improve the fit in the shorter wavelength we increase the WD temperature, 
as the mass accretion rate is already too large for quiescence and therefore cannot be increased. We 
find that the fit in the shorter wavelengths can be improved by elevating the WD to a temperature 
T ~ 50, OOOK, however the overall Xu is further increased by doing so. This model is shown in 
Figure 12. While the model does well at fitting the overall continuum of the spectrum, it does not 
fit the FUSE absorption lines as the velocity broadening in the disk makes the overall depth of the 
line shallower in the WD+disk model. It is possible that the absorption lines form above the disk 
and WD; we discuss this possibility and its implications in the next and last section. 



6. Discussion and Conclusions 



The over abundances and the unidentified lines in the FUSE spectrum of EM Cyg might be 
associated with some of the veiling material rather than the WD it self, we have no add itional way 

) show that 



to verify this. However, simulations of the stream disk interaction (jKunze et al 



2001 



the stream overflow hits the disk at smaller radii where it can in part bounce back. If this is the 
case, stream material can be found at any orbital phase, but with greater density and altitude 
above the disk plane around ~ 0.6 — 0.9. The rich variety of lines observed in the present FUSE 
spectrum of EM Cyg is reminiscent o f the HST/STIS s pectrum of the DN TT Crt observed in a low 
state with an inclination of z ~ 60° (ISion et al. II200S ), and of the FUSE spectrum of the NL BB 



Dor possibly with an inclination of 80° (jGodon et al. 



20081 ) . These spectra suggest line formation 



in different temperature regions, whether due to a wind above the disk (as could be the case in BB 
Dor) or overflowing stream material (possibly for TT Crt). So the strong and deep absorption lines 
in the FUSE exposure 1 of EM Cyg are likely due to circumstellar material, i.e. stream material 
overflowing the disk rim and bouncing off the disk at smaller radii, while the lower flux in exposures 
2 and 3 are due to the direct veiling of the WD by the stream overflow around phase 0.7. 



If the veiling material is responsible for the absorption lines, then the abundance fitting is 
not valid anymore, as the lines might be out of thermodynamic equilibrium or represent highly 
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variable optical depths. As a consequence, the best fit to the continuum is actually the WD+disk 
model shown in Figure 12, rather than the single WD model. While the single WD model has a 
lower value, the WD+disk model is in better agreement with the distance of the system and its 
higher-than-quiescent state. Because of this, we conclude that the FUSE spectrum of EM Cyg is 
best interpreted as the spectrum of a hot WD with a small contribution from a low mass accretion 
disk. 

The picture that emerges from this spectral analysis, is that of a 50,000K WD, contributing 
93% of the total FUV flux, and a disk, with M = 10~^'^MQyr~^ contributing the remaining 7%. 
Such a composite WD+disk model gives a distance d=382pc in best agreement with its expected 
value. The WD+disk model is also in agreement with the moderately higher-than-quiescent state 
during which the F{7S'£' spectrum was taken, even though the value is slightly larger as the model 
does not fit the absorption lines very well. The spectrum presents some deep absorption lines as 
well as unidentified lines in the longer wavelength region (> llOOA), apparently forming in front 
and above the WD/disk, from stream material overflowing the disk edge. As a consequence, the 
exposure taken at phase <j) ~ 0.45 has a higher flux than exposures obtained around (j) ~ 0.65 — 0.90. 
The large disk (0.3a) in EM Cyg is consistent with this picture, as the hot spot is consequently 
located closer to the secondary (around (j) ~ 0.9) and LI stream material can more easily be found 
off the disk plane. The large disk is also consistent with the non-negligible mass accretion rate 
(lO-iOMoyr-i). 

However, whether modeled with a single WD component or a WD+disk composite model, the 
FUSE spectrum of EM Cyg r eveals a rather hot WD- This result is i n stron g disagreement w i th the 



previous spectral analysis of IWelsh et al. I (|2005l ): IWinter Sion I (|2003l ): lUrban fc Sion I (|2006l ) 



who found a WD temperature T ~ 20 — 24, OOOK. The discrepancy with previous analysis is due 
to (i) veiling of the WD by stream material overflowing the disk rim (in both the FUSE and lUE 
spectra around orbital phase (p ~ 0.6 — 0.9, which were modeled as such); (ii) contribution of the 
hot spot and secondary impact region to the longer wavelengths of the lUE SWP spectra; and (iii) 
eclipse of the disk in the SWP 08088 lUE spectrum (taken around orbital phase 4> ~ 0.98 — 0.16). 
In addition, in these previous analysis, the flux was contributed mostly from an accretion disk, 
while in the present work the main component is the WD. 

Originally, the derivation of EM Cyg system parameters were affected by the presence of the 
third star, Arkadash, in the background of the system. In a similar way, we find that the previously 
derived mass accretion rate and WD temperature might have been affected by the veiling of the 
WD and disk by stream overflow around orbital phase (p ~ 0.65 — 0.90. 



PG wishes to thank Steve (Stephen) Lubow for some interesting discussions on the stream 
disk interaction, and Mario Livio for his kind hospitality at the Space Telescope Science Institute, 
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Association of Variable Star Observers (AAVSO) for providing public online optical archival data 
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Table 1. System Parameters for EM Cyg 



Parameter Value References 



Porb 7hr/0.29090912d Mumford k Krezminski ^1969) 

i 67 ± 2deg North et al. (2000) 

M^d IMq Welsh et al. (2007); Eobinson (1974); Stover et al. (1981) 

M2nd O.77M0 Welsh et al. (2007) ; Eobinson (1974) ; Stove r et al. (1981) 

d 350-500pc Bailev (1981); Jameson et al. (1981); Welsh et al. (2005) 

Vmin 14.4 Spodi et al. (2003, 2005) 

Knax 12.5 Spodi et al. (2003 . 2005 ^1 

DNO 14.6-16. 5s Stiening et al. (1982) 

2nd Spec.Type K3V North et al. (2000) 

Arkadash Spec.Type K2V-K5V North et al. (2000) 

Gal.Lat. b 4.277° SIMBAD 

E(B-V) 0.05±0.03 Verbunt ^987 ^1: La Dons (1991); Bruch fc Engel (1994) 
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Table 2. FUSE and lUE Observations Log 



Obs. 


Date 


Time 


Exp. time 


Dataset 


Start 


Phase^ 


V 




(dd/mm/yy) 


(UT) 


(sec) 




HJD 




(range) 


FUSE 


05/09/02 


llh34mlls 


1,479 


COlOOlOlOOl 


2452522.9849 


0.42-0.48 


13.1-13.9 


FUSE 


05/09/02 


12h54m33s 


2,651 


C0100101002 


2452523.0407 


0.62-0.73 


13.1-13.9 


FUSE 


05/09/02 


14h34m31s 


2,648 


C0100101003 


2452523.1102 


0.85-0.96 


13.1-13.9 


FUSE 


05/09/02 


16h22m59s 


2,135 


C0100101004 


2452523.1855 


0.11-0.19 


13.1-13.9 


lUE 


29/02/80 


14h42m50s 


2,700 


LWR 07056 


2444299.1102 


0.86-0.98 


13.6-13.7 


lUE 


29/02/80 


15h32m22s 


4,500 


SWP 08088 


2444299.1445 


0.98-0.16 


13.6-13.7 


lUE 


05/08/82 


02h40m47s 


1,500 


LWR 13862 


2445186.6152 


0.65-0.71 


13.2-13.5 


lUE 


05/08/82 


03hllm30s 


1,500 


SWP 17591 


2445186.6366 


0.72-0.78 


13.2-13.5 


lUE 


05/08/82 


03h44m26s 


2,400 


LWR 13863 


2445186.6595 


0.81-0.91 


13.2-13.5 


lUE 


05/08/82 


04h34m08s 


2,399 


SWP 17592 


2445186.6940 


0.92-0.02 


13.2-13.5 



^The ephemeris, time of minima /ecli pse HJD2400000 -|- 43780.7508, 45257.409L 
50692.4613, fc 53 7 09.18 73 were taken from iMnmfordI (|l98nl ): beuermann fc Paknll I (|l984l ): 
Csizmadia et al. I (|2008l ) with the period 0.29090912d=25134.55s 
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Table 3. WD Synthetic Spectra 





Vrotsini 


[Si] 


[S] 


[N] 


M 


i 


d 




WD/Disk 


Eb-v 


Figure 


(lO^K) 


(km/s) 


(Solar) 


(Solar) 


(Solar) 


(Moyr-i) 


(deg) 


(pc) 




% 










1.0 


1.0 


1.0 


10-9 


75 


441 


0.6916 


0/100 


0.00 


9 
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1.0 


1.0 


10-9 


60 
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0.6867 


0/100 


0.00 




35 


200 


1.0 


1.0 


1.0 






223 


0.6063 


100/0 


0.00 




40 


200 


1.0 


1.0 


1.0 






280 


0.5823 


100/0 


0.00 




40 
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1.0 


1.0 


1.0 
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0.5709 


100/0 


0.00 




40 


100 


30. 


10. 


1.0 
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0.4780 
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0.00 




40 
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10. 


10. 
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0.00 




40 
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264 


0.4768 


100/0 
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40 
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30. 
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10. 






261 


0.4767 


100/0 


0.00 


10 


40 
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30. 


10. 
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0.4776 
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0.05 




40 
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30. 


10. 


10. 
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60 
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0.5291 


77/23 


0.00 


11 
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60 
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Fig. 1. — The binary configuration of EM Cyg is shown during echpse of the disk at orbital phase 
= 0. The inchnation is i = 67° and the mass ratio is 0.77. The disk is echpsed but the white 
dwarf is not. 
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AAVSO DATA FOR EM CYG - WWW AAVSO.OHG 
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Fig. 2. — The visual magnitude of EM Cyg recorded by the members of the AAVSO 
(www.aavso.org). EM Cyg was observed with STIS (around JD2452515) when the system was 
in outburst (center of the graph), and was observed with FUSE a week later (on JD2452523) when 
the system was in quiescence. 
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Fig. 3. — The FUSE spectrum of EM Cyg. This spectrum consists of 4 exposures (correspond- 
ing to 4 FUSE orbits). The sharp emission hnes are contamination. The ISM molecular features 
are marked vertically for clarity. The molecular lines are identified by their band (W= Werner; 
L=Lyman), upper vibrational level (1-16), and rotational transition (R, P, or Q with lower rota- 
tional state J=l-3). The ISM atomic lines are marked below the x-axis. Some ISM metal lines 
are also identified: Ni, An. The Nitrogen lines are possibly also contaminated with terrestrial 
atmospheric absorption (especially the Ni IISSA). We have marked the position of the Niv 924 A 
lines, though they are not detected. The C iii and O vi lines are in broad emission, with a possible 
P-Cygni profile at the Cm 1175A . 
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Fig. 4. — The FUSE exposures 1 & 2 of EM Cyg. During the first exposure the flux is about 2 to 
4 times larger than during exposure 2. Exposure 3 is similar to exposure 2; exposure 4 is similar 
to exposure 1. For clarity exposures 3 and 4 are not shown, and the spectra have been binned at 
0.5A. 
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Fig. 5. — EM Cyg Binary Configuration during the 4 FUSE orbits/exposures. Exp.l (upper left) 
is at phase (j) ~ 0.45: the disk is in front of the observer and absorption is minimal. Exp. 2 (upper 
right) is at phase 4> ~ 0.65: the disk is still slightly in front moving away from the observer, the 
stream overflow is in front, absorption is maximal. Exp. 3 (lower left) is at phase (f) ~ 0.90: the 
secondary filling its Roche lobe is in front moving to the right, the stream and hot spot arc directly 
in front, absorption is maximal. Exp. 4 (lower right) is at phase (p ^ 0.15: the secondary Roche 
lobe is slightly in front and to the right of the observer, the disk is in the back moving towards the 
observer. The stream is in the back and behind, absorption is minimal. 




Fig. 6. — The Ff/S'E' exposures are shown together with the SWP /[/£' spectra. The FC/5'£' exposures 
1 and 4 have been combined together and are shown in red; exposures 2 and 3 have also been 
combined together and are shown in black. FUSE exposures 2 and 3 have a lower flux due to 
veiling of the WD by stream material overflowing the disk. It is likely that the lower flux in the 
lUE spectrum SWP 08088 is also due to the veiling of the WD. Additional change in the lUE 
spectra might be due to the cooling of the WD or ongoing accretion. 




Fig. 7. — The lUE LWR and SWP spectra of EM Cyg. The wavelength is in A, and the flux is 
in ergs s~^cm^^A^^. The plate au betwe e n 225 and 3000 A is possibly due to the hot spot with 
T ~ 10, OOOK, as pointed out bv lSzkodvl (|l98lh . 
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Fig. 8. — EM Cyg Binary Configuration during the 4 FC/S'i? orbits/exposures. The areas delimited 
in green represent the viewing angles during the 4 FUSE exposures (1,2,3 & 4). The disk radius 
here is 0.3a, where a is the binary separation. The stream trajectory from the LI point is shown 
in red. The stream overflows the disk rim and from there it is launched at an angle ~ 55° above 
the disk; it falls back onto the disk around (p ~ 0.5. The region for w hich the s trearn can reach 
z/r ~ 0.4 (corresponding to veiling at i ^ 67deg) is shown in blue (from Icodon I ( 2009 )). 
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Fig. 9. — The best fit single disk model (in solid black) to the FUSE/exposme 1 spectrum of EM 
Cyg. The observed spectrum (in red) has not been dereddened. The regions that have been masked 
for the fitting are shown in blue. The model consist of an accretion disk around a IMq WD with 
M = 1 X lO~^M0yr~-^, an inclination i = 75°. The distance obtained is d = 441pc and Xu ~ 0.6916. 
Because of the rotational broadening the model cannot fit the absorption lines from the source, 
and the mass accretion rate is too large and inconsistent with the quiescent state of the system. 
An ISM model has been added to the disk model to fit the ISM lines. The dotted line show the 
disk model without the inclusion of the ISM model. 
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Fig. 10. — The best fit WD model (in solid black) to the FUSE/ exposure 1 spectrum of EM Cyg. 
The observed spectrum (in red) has not been dereddened. The model consist of a I.OMq WD, with 
T = 40, OOOK, a projected rotational velocity Vrotsini = lOOkm/s, and over-solar abundances of Si 
(30), S (10), and N (10). The entire WD synthetic spectrum has been shifted to the blue by ~ O.SA, 
which corresponds to a radia l velocity of about ^ 150km/s (at 1,000A) in rough agreement with 
the radial velocity + 1 25km /s I Welsh et al. I (|2005l ). An ISM model (no shifted) has been added to 
improve the overall fit. The dotted line represents the WD model without the inclusion of the ISM 
model. The regions of the observed spectrum that have been masked for the fitting are shown in 
blue. The distance obtained from the best fit is d = 261pc and Xu — 0.4768 
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Fig. 11. — The best WD model (from Figm'e 10) plus a disk model are fitted (in solid black) to the 
FUSE/ exposure 1 spectrum of EM Cyg. The observed spectrum (in red) has not been dereddened. 
The regions that have been masked for the fitting are shown in blue. The disk model (dashed line) 
consist of an accretion disk around a IMq WD with M = 1 x lO~^^M0yr~^, an inclination i = 60°. 
The WD model (dotted line) consists of a l.OM© WD, with T = 40,000K, a projected rotational 
velocity K-ot sini = lOOkm/s, and over-solar abundances of Si (30), S (10), and N (10). The entire 
WD+disk synthetic spectrum has been shifted to the blue by ~ 0.5A . The distance obtained is 
d = 318pc and xl = 0.5291. An ISM model has been added to the WD+disk model to fit the ISM 
lines. The WD contributes 77% of the FUV flux and the disk contributes the remaining 23%. 
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Fig. 12. — Same as in Figure 11, but here the WD temperature has been increased to 50,000K. The 
WD contributes to 93% of the FUV flux and the disk contributes the remaining 7%. The distance 
obtained is d = 382pc and Xu — 0.5573. This model fits the shorter wavelength continuum flux 
better than the previous models, but it does not fit the absorption features as well as the single 
WD model. 



